RADIOWAVE MONITORING METHOD AND APPARATUS 

BACKGROUND OF THE INVENTION 

The present invention relates to a radiowave 
monitoring method and apparatus which can identify a 
location of a radiowave emitting source, such as an 
illegal radio station or the like, in a prescribed area 
and monitor radiowave using environments. 

Conventionally, in identifying a location of a 
radiowave emitting source, such as an illegal radio 
station in a prescribed area, monitor stations are 
disposed at a plurality of positions, and the respective 
radiowave stations observe arrival bearings of radiowaves 
from the radiowave emitting source by means of Yagi-Uda 
array antennas, goniometers or others. The respective 
monitor stations plot the arrival bearings on a map to 
thereby estimate the location of the radiowave emitting 
source, based on an intersection of the respective 
bearings . 

The above-described conventional method is based on 
the presumption that radiowaves from a radiowave emitting 
source being monitored propagate rectilinearly along a 
single path to a monitor station. Actually, however, 
radiowaves from a radiowave emitting source propagate to 
a monitor station, in some cases being diffracted by or 
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reflected on topographies, and ground objects, such as 
buildings, etc. Radiowaves from a radiowave emitting 
source propagate to a monitor station, in other cases, in 
multi-paths under influence of topographies and ground 
objects. The estimation of a location of a radiowave 
emitting source being monitored by the conventional 
method has a problem that radiowaves propagate, being 
diffracted or reflected or in multi-paths, which lowers 
accuracy of estimating the location of the radiowave 
emitting source. 

The inventors of the present application already 
proposed a wide-area radiowave monitoring method and 
apparatus which can solve the above-described problem 
(refer to, e.g., Japanese Patent Application Laid-Open 
Publication No. Hei 11-326482 (1999)). 

In the wide-area radiowave monitoring method 
described in Japanese Patent Application Laid-Open 
Publication No. Hei 11-326482 (1999), in a computer 
simulation, one or more monitor stations are disposed in 
a certain area, radiowave emitting sources are present at 
a plurality of locations in the area, and arrival 
bearings of radiowaves from the radiowave emitting 
sources, which income to the monitor stations are 
computed for the respective radiowave emitting sources, 
using topographical information, while the monitor 
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stations observe for respective divided propagation paths 
radiowaves from a radiowave emitting source. The arrival 
bearings observed by the monitor stations, and results of 
the computer simulation are compared with each other to 
probe the computer simulation results for arrival 
bearings of the results, which are most similar to the 
observed arrival bearings of the radiowaves, and a 
corresponding location is judged to be an alleged 
location of the radiowave emitting source. Further, on 
the presumption that the radiowave emitting source is 
present near the alleged location, the computer 
simulation is executed, adjusting a location of the 
radiowave emitting source, to thereby determine a 
location of the radiowave emitting source. 

The above-described method can accurately determine 
a location of a radiowave emitting source even in a case 
that radiowaves are diffracted or reflected. 

However, in the wide-area radiowave monitoring 
method and apparatus described in the specification of 
Japanese Patent Application Laid-Open Publication No. 
326482/1999, arrival bearings of radiowaves incoming to 
the monitor stations from radiowave emitting sources at a 
plurality of locations in an area are computed by the 
computer simulation for the respective locations, which 
requires a plurality of sets of a data base provided by 
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the computer simulation of the radiowave propagation to 
be prepared for unidentified radiowave emitting sources. 
Furthermore, it takes much time to prepare the data base. 



SUMMARY OF THE INVENTION 
An object of the present invention is to provide a 
radiowave monitoring method and apparatus which can 
shorten time for preparing a data base for the radiowave 
propagation simulation to make the radiowave monitor more 
efficient. 

The above-described object is achieved by a 
radiowave monitoring method comprising: the step of 
comparing an observation result of an arrival direction 
and a pattern of intensities of a radiowave observed at 
one position in an observation area with simulation 
results of arrival directions and patterns of intensities 
of radiowaves emitted from said one position, which are 
to be observed at other plural positions in the 
observation area. 

In the above-described radiowave monitoring method, 
it is possible that of said plural positions whose 
simulation result shows the arrival direction and the 
pattern of intensities of the radiowave most correlated 
with the observation result at said one position is 
identified as a location of a radiowave emitting source. 

- 4 - 



In the above-described radiowave monitoring method, 
it is possible that in simulating the arrival directions 
and the patterns of intensities of radiowaves emitted 
from said one position, the observation area is 
two-dimensionally divided into a plurality of regions, 
and electric field intensities to be observed in the 
respective regions are computed. 

In the above-described radiowave monitoring method, 
it is possible that in simulating the arrival directions 
and the patterns of intensities of radiowaves emitted 
from said one position, the observation area is 
three-dimensionally divided into a plurality of spaces, 
and electric field intensities to be observed in the 
respective spaces are computed. 

In the above-described radiowave monitoring method, 
it is possible that in simulating the arrival directions 
and the patterns of intensities of radiowaves emitted 
from said one position, changing the emission direction 
of radiowaves emitted from said one position, the 
electric field intensities to be observed in the 
respective regions or the respective spaces are computed 
to give electric field intensity distributions for the 
respective emission directions. 

In the above-described radiowave monitoring method, 
it is possible that in simulating the arrival directions 
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and the patterns of intensities of radiowaves emitted 
from said one position, geography and ground objects in 
the observation area are taken into consideration. 

In the above-described radiowave monitoring method, 
it is possible that based on the location of the 
identified radiowave emitting source and the simulation 
result, propagation path of the radiowaves from said one 
position to the radiowave emitting source is traced. 

In the above-described radiowave monitoring method, 
it is possible that based on a result of tracing the 
propagation path, antenna directivity of the radiowave 
emitting source is estimated. 

In the above-described radiowave monitoring method, 
it is possible that based on the estimated antenna 
directivity of the radiowave emitting source, an electric 
field intensity distribution of the radiowaves emitted 
from the radiowave emitting source is computed. 

The above-described object is also achieved by a 
radiowave monitoring apparatus comprising: a radiowave 
observing means disposed at one position in an 
observation area, for observing arrival directions and 
patterns of intensities of radiowaves; a storing means 
for storing simulation results of arrival directions and 
patterns of intensities of radiowaves emitted from said 
one position, which are to be observed at other plural 
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positions in the observation area; and a radiowave 
emitting source identifying means for comparing an 
observation result given by the radiowave observing means 
with the simulation results stored in the storing means 
to identify that of said plural positions whose 
simulation result shows the arrival direction and the 
pattern of intensities of the radiowave most correlated 
with the observation result at said one position as a 
location of a radiowave emitting source. 

In the above-described radiowave monitoring 
apparatus, it is possible that the storing means 
two-dimensionally divides the observation area into a 
plurality of regions, and computes electric field 
intensities of radiowaves emitted from said one position, 
which are to be observed in the respective regions. 

In the above-described radiowave monitoring 
apparatus, it is possible that the storing means 
three-dimensionally divides the observation area into a 
plurality of spaces and computes electric field 
intensities of radiowaves emitted from said one position, 
which are to be observed in the respective spaces. 

In the above-described radiowave monitoring 
apparatus, it is possible that the storing means stores 
electric field intensities of radiowaves emitted from 
said one position in different directions, which are to 
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be observed in the respective regions or the respective 
spaces, for the respective directions. 

In the above-described radiowave monitoring 
apparatus, it is possible that the apparatus further 
comprises a propagation path tracing means for tracing a 
propagation path of radiowaves from said one position to 
the radiowave emitting source, based on the location of 
the radiowave emitting source identified by the radiowave 
emitting source identifying means and the simulation 
result. 

In the above-described radiowave monitoring 
apparatus, it is possible that the apparatus further 
comprises an antenna directivity estimating means for 
estimating antenna directivity of the radiowave emitting 
source, based on a result of tracing the propagation path 
given by the propagation path tracing means. 

In the above-described radiowave monitoring 
apparatus, it is possible that the apparatus further 
comprises an electric field intensity computing means for 
computing an electric field intensity distribution of 
radiowaves emitted from the radiowave emitting source, 
based on the antenna directivity of the radiowave 
emitting source estimated by the antenna directivity 
estimating means . 

As described above, according to the present 
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invention, an observation result at one position is 
compared with simulation results of arrival directions 
and patterns of intensities of radiowaves emitted from 
said one position, which are to be observed at other 
plural positions, and that of said plural positions whose 
simulation result shows the arrival direction and the 
pattern of the intensities of the radiowave most 
correlated with the observation result at said one 
position is identified as a location of the radiowave 
emitting source. Accordingly, a time for preparing data 
base by the radiowave propagation simulation can be 
decreased, and the radiowave monitor can be more 
efficient . 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagrammatic view showing the radiowave 
monitoring method according to one embodiment of the 
present invention. 

FIG. 2 is a block diagram of a constitution of the 
radiowave monitoring apparatus according to the 
embodiment of the present invention. 

FIG. 3 is a diagrammatic view of a radiowave 
propagation simulation model. 

FIG. 4 is a diagrammatic view of a coordinate system 
for computing diffraction coefficients. 
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FIGs. 5A-5G are views explaining an estimation of 
directivity of a transmission antenna of a radiowave 
emitting source. 

FIG. 6 is a flow chart of a process of the radiowave 
monitoring method according to the embodiment of the 
present invention. 

FIG. 7 is a diagrammatic view of a three-dimensional 
radiowave propagation simulation model. 

DETAILED DESCRIPTION OF THE INVENTION 
The radiowave monitoring method and apparatus 
according to one embodiment of the present invention will 
be explained with reference to FIGs. 1 to 6 . FIG. 1 is a 
diagrammatic view of a constitution of the radiowave 
monitoring method according to the present embodiment. 
FIG. 2 is a block diagram of a constitution of the 
radiowave monitoring apparatus according to the present 
embodiment. FIG. 3 is a diagrammatic view of a radiowave 
propagation simulation model. FIG. 4 is a diagrammatic 
view of a coordinate system for computing diffraction 
coefficients. FIGs. 5A-5G are views explaining an 
estimation of directivity of a transmission antenna of a 
radiowave emitting source. FIG. 6 is a flow chart of the 
process of the radiowave monitoring method. 
(Radiowave Monitoring Apparatus) 
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The radiowave monitoring apparatus according to the 
present embodiment will be summarized with reference to 
FIG. 1. 

A radiowave emitting source 10 is an emitting source 
for emitting radiowaves to be monitored. Here the 
radiowave emitting source 10 emits radiowaves of a 
frequency f 1 . For radiowave monitor of radiowaves emitted 
by the radiowave emitting source 10, monitor stations 12a, 
12b, 12c which make holographic observation of emitted 
radiowaves are positioned in a monitor area. The sensor 
stations 12a, 12b, 12c are connected to a center station 
16 via circuits 14, which performs radiowave propagation 
simulation and compares simulation results with 
observation results to identify a location of the 
radiowave emitting source 10. 

Each monitor station 12a, 12b, 12c has a radiowave 
hologram observation and reconstruction apparatus 18 for 
making holographic observation of radiowaves and 
outputting reconstructed radiowave images. The radiowave 
hologram observation and reconstruction apparatus 18 
includes, e.g., a fixed antenna, and a scanning antenna 
which scans a scan observation plane. Radiowave images 
are reconstructed by correlating for a prescribed 
observation frequency received signals received by the 
fixed antenna and received signals received by the 
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scanning antenna with each other to give correlated 
values (two-dimension complex inter-holograms) and 
reconstructing the two-dimensional inter-holograms. 

Radiowaves from the radiowave emitting source 10 
propagate to the respective monitor stations 12a, 12b, 
12c. The radiowaves contain, in addition to a component 
(direct waves) which have propagated from the radiowave 
emitting source 10 directly to the monitor stations, a 
component (diffracted waves) and a component (reflected 
waves) which have propagated, respectively diffracted and 
reflected by mountains and other ground objects. In the 
present embodiment, each monitor station 12a, 12b, 12c 
makes radiowave hologram observation and reconstructs 
radiowave images so as to separate components of 
respective propagation paths to give their arrival 
directions and intensities. Radiowaves, which are one 
kind of waves, can be holographically observed in the 
same way as in the holographic observation of light; 
holograms are reconstructed to obtain reconstructed 
radiowave images, and based on the reconstructed 
radiowave images, a wave source distribution, intensities, 
etc. can be investigated. When component of respective 
propagation paths can be separated to be observed, the 
conventional direction finding technique may be used in 
place of the radiowave hologram observation. 
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The center station 16 includes simulation means 20 
for performing computer simulation of radiowave 
propagation, a data base 2 6 for storing a result of a 
simulation made in advance by the simulation means 20, 
and a comparing unit 28 which compares reconstructed 
radiowave images observed by the respective monitor 
stations 12a, 12b, 12c with an electric field intensity 
distribution given by the simulation. The center station 
16 further includes a direction finding trace drawing 
preparing unit 3 0 which identifies a location of the 
radiowave emitting source 10, based on a comparison 
result given by the comparing unit 28 and prepares to 
output a direction finding trace drawing, and a radiowave 
activity map preparing unit 32 which prepares and outputs 
in the form of a radiowave activity map an electric field 
intensity distribution of radiowaves from the radiowave 
emitting source 10 whose location has been identified. 

The simulation means 20 includes a map information 
storing unit 22 which stores information (map 
information) of a map, geography, ground objects, etc., 
and a radiowave propagation simulator 24 which simulates 
with reference to the map information radiowave emission 
from the monitor station 12 on the map to simulate 
radiowave propagation. The radiowave propagation 

simulator 24 makes a computer simulation of an electric 
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field intensity distribution generated when radiowave 
emission from a monitor station 12 on the map data is 
simulated, by using, e.g., ray tracing (refer to, e.g., 
IEEE Network Magazine, pp. 27-30, November, 1991) or 
moment method (refer to, e.g. r.f. Harrington, Field 
Computation by Moment Methods, IEEE Press, 1993). 

The data base 28 stores a matrix of the electric 
field intensity distribution computed by the simulation 
means 20. The matrix of the electric filed intensity 
distribution stored by the data base 26 is referred to by 
the comparing unit 2 8 when a location of the radiowave 
emitting source 10 is identified. 

The comparing unit 2 8 compares a computer simulation 
result of the electric field intensity distribution for 
the simulated radiowave emission from the monitor station 
12a, which have been computed by the simulation means 20 
and are stored by the data base 26 with a result of 
radiowave holographical observation made by the actual 
monitor station 12a. Based on the result of the 
comparison, a location of the radiowave emitting source 
10 is identified. 

The direction finding trace drawing preparing unit 
30 traces radiowave propagation paths from the radiowave 
emitting source 10 whose location has been identified by 
the comparing unit 28 to the respective monitor stations 
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12a, 12b, 12c to prepare and output a direction finding 
trace map 34, and furthermore estimates an antenna 
directivity of the radiowave emitting source 10, based on 
result of the trace. 

Based on the location of the radiowave emitting 
source 10 identified by the comparing unit 28, and an 
estimation result of the antenna directivity of the 
radiowave emitting source 10 given by the direction 
finding trace drawing preparing unit 30, the radiowave 
activity map preparing unit 32 commands the simulation 
means 20 to execute the computer simulation of radiowave 
propagation to compute an electric field intensity 
distribution of radiowaves from the radiowave emitting 
source 10 on the map. A computed result is displayed and 
printed out in a radiowave activity map 3 6 of the 
radiowave emitting source 10. 

The radiowave monitoring apparatus according to the 
present embodiment, which has the above-described 
constitution, is characterized in that for simulated 
radiowave emission from the monitor station 12a on the 
map an electric field distribution is computed in advance 
by computer simulation for respective regions 
corresponding to direction of the radiowave emission. 
Arrival directions and intensities of radiowaves given by 
reconstructed images of actual radiowaves observed by the 
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monitor station 12a are compared with the electric field 
intensity distribution given by the computer simulation 
to thereby identify a location of the radiowave emitting 
source 10. 

(Radiowave Monitoring Method) 

Next, the radiowave monitoring method according to 
the present embodiment will be explained with reference 
to FIGs. 1 to 6. 

Before electromagnetic waves are actually monitored, 
for simulated radiowave emission from the monitor station 
12a on map data an electric field intensity distribution 
is computer simulated in two-dimensional model. The 
computer simulation for the electric field intensity 
distribution will be detained with reference to FIGs. 3 
and 4. The computer simulation which will be detailed 
below uses ray tracing. The ray tracing is a method in 
which it is assumed that a number of rays (light rays) 
are emitted by a radiowave emitting source, and the rays 
are traced along their paths at certain angles. When a 
ray hits the ground or a building, the ray is further 
traced in a reflected direction from the position as a 
reflection point. When a ray hits an edge of a building, 
a plurality of rays are generated at the point as a 
diffraction point. 

First, the radiowave propagation simulator 24, with 
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reference to map data of an area for the simulation to be 
executed read from the map information storing unit 22, 
divides the area for the simulation to be executed into 
an M x N matrix of regions. The respective divided 
regions are represented by 1, 2, . .., MN-1, mn. 
Information of geography and ground objects is mapped in 
the map data- In FIG. 3, a building 13a and a building 
13b are present. 

Subsequently, as shown in FIG. 3, a ray is emitted 
in an L k direction from the monitor station 12a, and 
propagation paths of the ray and electric field 
intensities of the ray passing through the respective 
regions are computed. Computed electric field intensities 
of the ray are recorded as elements of the matrix having 
ray numbers L k of emitted rays as the row numbers and 
region numbers as the column numbers . The matrix for 
recording the electric field intensities is shown in 
TABLE 1 below. 



TABLE 1 




1 



2 



3 



MN 



1 



En E 3 



'12 



E 3 



2 




3 



At this time, the ray has a thickness proportional 
to a propagation distance from the monitor station 12a as 
an origin, and when the ray of the thickness covers parts 
of the respective regions, the ray is considered to have 
passed the regions. When a plurality of rays have passed 
through one region, an electric field intensity of that 
of the rays which gives a highest intensity is recorded 
in the matrix of TABLE 1 . 

Such simulation process will be explained 
specifically with reference to a case as shown in FIG. 3 
in which a ray emitted from the monitor station 12a as an 
origin in a direction L k is reflected on the building 13a 
at a point p and arrives at the building 13b at a point q. 
Here, a line segment interconnecting the monitor station 
12a and the point p is represented by a line segment a, 
and the line segment interconnecting the point p and the 
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point q is represented by a line segment b. 

For regions through which the line segment a of a 
ray emitted from the monitor station 12a as an origin to 
the point p has passed, electric intensities computed by 
the following Formula 1 are recorded in corresponding 
elements <L k/ region numbers) of the matrix. 

*3S=M 

r x 1 

wherein E 0 represents an electric field of radiowaves when 
emitted by the monitor station 12a, r represents a 
distance from the monitor station 12a to the centers of 
the respective regions, and k is a wave number of 
radiowaves . 

Subsequently, for regions through which the line 

segment b which has been reflected on the building 13a at 

the point p and arrived at the point q of the building 

13b, electric field intensities computed by the following 

Formula 2 are recorded in corresponding elements (L k , 

region numbers) of the matrix. 

exp(-jka) a 

E 0 ^ J -R-exp(-)kr') (2 ) 

r a + r' 

wherein a represents a length of the line segment a, and 

r" represents a distance from the point p as an origin to 

the centers of the respective regions. 

R represents a reflection coefficient and is given 

as follows. 
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First, a complex index of refraction n at the time 
when a plane wave is incident on a medium 2 from a medium 
1 in accordance with Snell's law is expressed by the 
following Formula 3 . 



(3) 



wherein z ir ^ and o 1 respectively represent a dielectric 
constant, a magnetic permeability and a conductivity of 
the medium 1. e 2 , \i 2 and o 2 respectively represent a 
dielectric constant, a magnetic permeability and a 
conductivity of the medium 2. A reflection coefficient R 
is expressed by the following Formula 4 using a complex 
index of diffraction n when an incident electric field is 
in incidence plane, i.e., in a case of TM (Transverse 
Magnetic) incidence. 



R ^n 2 cosfl - [A 2 jn 2 - sin 2 8 



2 n — —r- ( 4 > 

/UjH cos# + ,u 2 Vh -sin 6 
On the other hand, when an incident electric field 
is vertical to incidence plane, i.e., in a case of te 
(Transverse Electric) incidence, a reflection coefficient 
R is expressed by the following Formula 5. 



R 



^ cosd - ii 2 4n z -sin 2 0 



cos# + fj, 2 ^n 2 - sin 2 0 ^ ^ 

Subsequently, line segments c 17 c 2 , c 3 , ... of 
diffracted waves of the ray which has been reflected on 
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the building 13a at the point p and arrived at the point 
q of the building 13b are given. Here, when the ray is 
incident on a corner of the building 13b, i.e., the point 
q is positioned at the corner of the building 13b, a 
plurality of rays are surface radiated from the corner. 

For regions through which the line segments c 2 , c 2 , 

c 3 , have passed, the following Formula 6 is recorded 

in corresponding elements (L k , region numbers) of the 
matrix. 



a+b Fy J } \r"(r"+a+b) ^ } 

wherein b represents a length of the line segment b, r" 
represents distances from the point q as an origin to the 
centers of the respective regions, and D represents a 
diffraction coefficient, and is computed by the following 
Formula 7, using the coordinate system shown in FIG. 4. 

D = -f}< 0 P 0 D s -j>'i>D h (7 ) 
wherein p 0 , p 0 ' , <j>, (j>', p d and S d are physical parameters 
indicative of the positional relationships shown in FIG. 
4. In FIG. 4, E° represents an incident wave, and E d 
represents a diffracted wave. D s and D h are expressed by 
the following Formula 8 which uses UTD (refer to R.G. 
Kouyoumjian, P.H. Pthak, Proceedings of the IEEE, no. 11, 
pp. 1448-1461, December, 1974) using Fresnel integration 
function. 
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-exp 




\F(kLa + (<t>-<f><)) 



sin/? 0 \[ [ In j 




]F(kLa- (</>-$')) 




\F(kLa\4> + <t><)) 




Here, L, a + (x) and a~(x) are respectively expressed 



by the following Formula 9. 




(9) 



N + and N are expressed by the following Formula 10. 



The above-described operation is repeated, changing 
a direction L k of the ray emitted by the monitor station 
12a to complete the matrix of an electric field intensity 
distribution for L k shown in TABLE 1. The matrix of the 
thus-given electric field intensity distribution is 
recorded in the data base 26. Thus, the simulation of 
radiowave propagation is completed. 

In a case that, as shown in FIG. 1, a plurality of 
monitor stations 12a, 12b, 12c are positioned, the above- 
described simulation is performed by the monitor stations 




(10) 
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besides the monitor station 12a, whereby the monitor can 
be made simultaneously by a plurality of monitor stations 
A plurality of monitor stations are used, whereby no 
insensitive region is present, and in addition, increased 
monitored data can improve accuracy of estimating a 
location of a radiowave emitting source. 

As described above, the radiowave propagation 
simulation is made in advance before the radiowave 
monitor is started. The radiowave monitor is made in the 
sequence of the flow chart shown in FIG. 6. 

First, the radiowave hologram observation is 
performed by the monitor station 12a to obtain 
reconstructed radiowave images (Step Sll). The radiowave 
hologram observation enables the measurement to be made 
in correlation between arrival directions and intensities 
of radiowaves . 

Then, the comparing unit 2 8 compares arrival 
directions and intensity patterns of radiowaves based on 
the reconstructed radiowave images observed by the 
monitor station 12a with respective patterns of column 
vectors of the matrix of electric field intensity 
distributions stored in the data base 26 (Step S12). 

In radiowave propagation, generally a raiowave 
propagation path is reversible between the emission and 
the receipt, and reversibility of a propagation 
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attenuation amount is true, i.e., a propagation 
attenuation amount is the same reversibly in the emission 
and in the receipt. 

Accordingly, out of column vectors of the matrix of 
electric field intensity distributions, a region having a 
region number of a most identical column vector is a 
location of the radiowave emitting source 10. 

The use of the method for thus identifying the 
radiowave emitting source 10 makes it unnecessary that 
when it is assumed that a radiowave emitting source is 
located at one of a plurality of positions in an area to 
be radiowave-monitored, reconstructed radiowave images 
observed by a monitor station are computed for the 
respective points by the computer simulation using map 
information, which can improve efficiency of the 
radiowave monitor. 

For the above-described comparison for identity, the 
usual pattern matching may be used. For example, when a 
highest level of radiowaves of a reconstructed radiowave 
image is 0 dB, patterns for radiowave levels above -40 dB 
are compared to give a correlation coefficient as shapes. 
As an antenna directivity of a radiowave emitting source 
10 is unknown, higher identity in arrival directions of 
radiowaves may be used. 

When the monitor stations 12b and the monitor 
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stations 12c in addition to the monitor station 12a 
concurrently perform the monitor, a result of a simulated 
electric field intensity distribution of radiowave 
emission simulated in advance is stored in the data base 
26 to be compared with a result given by actual radiowave 
hologram observation. 

Next, after a location of the radiowave emitting 
source 10 is identified, radiowave propagation paths from 
the location of the radiowave emitting source 10 
identified by the comparing unit 2 8 to the respective 
monitor stations 12a, 12b, 12c are traced by the 
direction finding trace preparing unit 30 to prepare and 
output a direction finding trace drawing 34 (Step S13). 

Path information of rays corresponding to the matrix 
of the given electric field intensity distribution is 
recorded in the data base 26 when the computer simulation 
of an electric field intensity distribution of simulated 
radiowave emission from the monitor station 12a on the 
map data is executed, whereby simultaneously with the 
identification of a location of a radiowave emitting 
source 10, a result of tracing radiowave propagation 
paths can be given without the use of the direction 
finding trace drawing preparing unit 30. 

Furthermore, based on a result of the above- 
described trace, antenna directivity of the radiowave 
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emitting source 10 is estimated as follows (Step S14). 

First, each monitor station 12a, 12b, 12c perform 
radiowave hologram observation to obtain reconstructed 
radiowave images. Arrival directions and intensities of 
the radiowaves for respective propagation paths are 
extracted from the reconstructed radiowave images 
obtained by each monitor station 12a, 12b, 12c. 

FIGs. 5A to 5C respectively show amplitudes of the 
radiowaves observed by the monitor stations 12a, 12b, 12c. 
The thick arrows indicate arrival directions of 
radiowaves along propagation paths, and the lengths of 
the arrows indicate amplitudes of the radiowaves. In the 
shown example, amplitudes for the respective propagation 
paths are a ir a 2 , h ir b 2 , c lf c 2 . On the other hand, 
radiowave propagation is simulated by a simulator 23 on 
the assumption that a non-directional transmitting 
antenna is present at the location of the radiowave 
emitting source 10, which has been already identified, to 
give amplitudes of the radiowaves for the respective 
propagation paths by the respective monitor stations 12a, 
12b, 12c. 

FIGs. 5D to 5F show amplitudes of radiowaves 
simulated by the respective monitor stations 12a, 12b, 
12c. Amplitudes for respective propagation paths are a 1 ' , 
a 2 ' , b 1 ' , b 2 ' , c x ' , and c 2 ' . Amplitudes observed for the 
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respective propagation paths P lf P 2 , p 3 , p 4f p 5/ p 6 are 
divided by amplitudes given by the simulation to give 
directivity of the transmitting antenna along the 
respective propagation paths at the location of the 
radiowave emitting source 10. 

In FIG. 5G, the black points indicate directivity 
along the respective propagation paths. Interpolation 
computation is made in consideration of an estimated kind 
of the antenna to thereby compute directional 
characteristics for arbitrary directions. The broken line 
in FIG. 5G indicates the thus-computed directional 
characteristics . 

Based on a result of the estimation of an antenna 
directivity of the radiowave emitting source 10 by the 
direction finding trace drawing preparing unit 30, the 
radiowave acitivity map preparing unit 32 commands the 
simulation means 20 to execute the computer simulation of 
radiowave propagation to compute an electric field 
intensity distribution of the radiowaves from the 
radiowave emitting source 10 on a map. A result of the 
computation is displayed and printed out as a radiowave 
activity map 36 of the radiowave emitting source 10 (Step 
S15) . 

Thus, the radiowave emitting source 10 in a certain 
area is identified, and at the same time, a direction 
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finding trace drawing is prepared, a directivity of the 
transmitting antenna of the radiowave emitting source 10 
is estimated, and a radiowave activity map is prepared, 
and the radiowave monitor is completed. 

As described above, according to the present 
embodiment, when it is assumed that a radiowave emitting 
source is present at one of a plurality of positions in 
an area, it is not necessary to compute by the computer 
simulation reconstructed radiowave images observed by the 
monitor station, which can decrease time for preparing 
data base and accordingly can improve radiowave 
monitoring efficiency. 

Next, the radiowave monitoring method and apparatus 
according to another embodiment of the present invention 
will be explained with reference to FIG. 7. FIG. 7 is a 
diagrammatic view of a three-dimensional radiowave 
propagation simulation model. 

The radiowave monitoring method and apparatus 
according to the present embodiment is characterized in 
that an electric field intensity distribution of 
simulated radiowave emission by a monitor station 12a on 
map data is simulated by using a three-dimensional model. 
The present embodiment is the same as the above-described 
embodiment except that an electricl field intensity 
distribution is simulated by using a three-dimensional 
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model . 

First, as shown in FIG. 7 , a monitor area is divided 
in x x y x z unit spaces, and the respective spaces are 
represented by space numbers 1, 2, 3, xyz-1, xyz . 

Subsequently, one ray is emitted in an L k direction 
from a monitor station 12a on map data. Propagation paths 
of the ray and electric field intensities of the ray at 
the time when the ray passes the respective unit spaces 
are computed in the same way as in the two-dimensional 
model of the above-described embodiment. Computed 
electric field intensities are recorded as elements of 
the matrix of emitted ray numbers L k in the rows and unit 
space numbers in the columns . 

At this time, the ray has a thickness proportional 
to a propagation distance from the monitor station 12a as 
an origin, and when the ray of the thickness covers parts 
of the respective regions, the ray is considered to have 
passed the regions. When a plurality of rays have passed 
through one region, an electric field intensity of that 
of the rays which gives a highest intensity is recorded 
in the matrix. 

When a ray emitted from the monitor station 12a in 
an L k direction is reflected on a building 13a at a point 
p or diffracted at a point q of building 13b, electric 
field intensities of the ray are computed in the same way 
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as in the above-described embodiment. Directions of the 
reflection and the diffraction are determined in 
consideration of three-dimensional information of the 
buildings 13a, 13b. 

The above-described operation is repeated, three- 
dimensionally changing directions of rays emitted by the 
monitor station 12a to thereby complete the matrix of an 
electric field intensity distribution of rays L k in the 
same way as in the two-dimensional model. The thus- 
obtained matrix of an electric field intensity 
distribution is recorded in the data base 26. Thus, the 
radiowave propagation simulation is completed. 

Using a result of the above-described radiowave 
propagation simulation, the radiowave monitor can be 
performed in accordance with the flow chart shown in FIG. 
6 in the same was as in the above-described embodiment. 

As described above, according to the present 
embodiment, it is unnecessary that when it is assumed 
that a radiowave emitting source is located at one of a 
plurality of positions in an area to be radiowave- 
monitored, reconstructed radiowave images observed by a 
monitor station are computed for the respective points by 
the computer simulation using map information, which can 
decrease time of preparing data base and can improve 
efficiency of the radiowave monitor. 
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